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ABSTRACT 

We report on the first detection of ionospheric disturbances caused by short 
repeated gamma-ray bursts from the magnetar SGR J1550— 5418. Very low fre- 
quency (VLF) radio wave data obtained in South America clearly show sudden 
amplitude and phase changes at the corresponding times of eight SGR bursts. 
Maximum amplitude and phase changes of the VLF signals appear to be corre- 
lated with the gamma-ray fluence. On the other hand, VLF recovery timescales 
do not show any significant correlation with the fluence, possibly suggesting that 
the bursts' spectra are not similar to each other. In summary, the Earth's iono- 
sphere can be used as a very large gamma-ray detector and the VLF observations 
provide us with a new method to monitor high energy astrophysical phenomena 
without interruption such as Earth Occultation. 

Subject headings: neutron stars: general — neutron stars: individual(SGR J1550- 
5418) 
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Introduction 



Very low frequency (VLF; 3—30 kHz) radio waves are reflected at the Earth' s lower iono - 



sphere and ground, and propagate within the Earth-ionosphere waveguide (e.g.. lWaitlll962l ). 
Since amplitude and phase of VLF radio waves are sensitive to the condition of the lower iono- 
sphere, they have been utilized to investigate the physics of the lower ionosphere. Energetic 
electrons can precipitate into the ionosphere due to wave-particle interact ion in the magne- 
tosph ere and cause VLF signal amplitude and phase perturbations (e.g., iKikuchi fc Evans 
19831 ). Soft X-rays from solar flares are a lso another source o f ionospheric distur bances 



which are detected using VLF signals (e.g., iTodoroki et al.l 120071 ; iRaulin et al.ll2010l ). 



Besides these solar-terrestrial events, the lower ionosphere is also affected by high energy 
photons (X-rays and gamma-rays) from extra-solar sources. An ionospheric d isturb ance 
caused by a cosmic gamma-ray burst was first reported by iFishman fc Inanl (119881 ) . It 



suggested that gamma-rays deposit their energies in the lower ionosphere, ionize abnormally 
the neutral atmosphere there, and modify the electron density height profile. In addition, 



Inan et al. 


1999, 


2007; 


Tanaka et al. 


2008) 



Magnetars emit a lot of shor t-duration gamma-ray flares repeatedly during active phases 



Woods fc Thompson! 120061 ). Typical duration and flux of the short repeated bursts are 
0.1 — 1 s and 10 -6 erg s _1 cm -2 , respectively. Furthermore, magnetars rarely emit excep- 
tionally bright gamma-ray flares (giant flares). So far, only three g iant flares were rec orded. 
The first one was detected in 1979, from the source SGR 0526—66 (IMazets et al.lll979l ). The 
second and third ones were emitted by SGR 1900+14 and SGR 1806—20, and they wer e ob- 



serve d by satellites in 1998 (IHurlev et a 



2004 ( ITerasawa et al.ll2005l : IHurlev et al. 



1999 



2005; 



Mazets et a 



Palmer et al. 



1999 



2005 



Tanaka et al 



Frederiks et al. 



200711 and 



20071 1. re- 



spectively. Since the fluences were much larger than those of GOES X-class solar flares by 
a few orders of magnitude, ionospheric disturbances caused by the se giant flares were de- 
tected as sudden a nd large amplitude changes of VLF radio waves (jlnan et al.lll999l 120071 : 
Tanaka et al~]l2008h . 



On the other hand, VLF amplitude and phase changes caused by short repeated gamma- 
ray bursts from a magnetar have not been detected so far because of the lack of high sen- 
sitivity of VLF observing systems. On 22 January 2009, one of the know n magnetars SGR 



J155 0— 5418 emitted a lot of short-duration gamma-ray bursts repeatedly (IMereghetti et al. 



20091 ). In this paper, we report on the first VLF detection of short repeated gamma-ray flares 
from this object. VLF data were pro vided by the South America VLF Network (SAVNET) 
tracking system ( IRaulin et al.l 120091 ). In §2, we describe details of the SAVNET observa- 
tions. Comparison of VLF amplitude and phase changes with gamma-ray fluences measured 
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by INTEGRAL satellite are presented and discussed in §3 and §4. We summarize this paper 
in 55. 



2. Observations 

The VLF data shown in this paper were obtained by SA VNET, which was recently in- 



stalled in Brazil, Peru, and Argentina (see lRaulin et al.l (120091 ) for the details of the SAVNET 



instrumental facility). Figure 1 shows locations of a relevant observing station ATI (Atibaia, 
Sao Paulo, Brazil) as well as five VLF transmitters (NPM, NLK, NDK, NAU, and NAA), 
and VLF signals from them have been continuously recorded with the time resolution of 1 
sec. The propagation path from NPM (21.4 kHz) to ATI is also drawn in Figure 1. Shaded 
hemisphere in Figure 1 exhibits the night-side part of the Earth at 6:48 UT, when the most 
intense gamma-ray flare occurred. The point on the Earth directly beneath the flare (sub- 
flare point) was located at 54.3°S, 14.0°E, and its position is shown using a cross. The part 
of the Earth illuminated by gamma-rays at 6:48 UT is illustrated by dashed area. 

Figure 2 shows NPM-ATI amplitude and phase data recorded from 4:00 UT to 10:00 
UT on 22 January, 2009. In Figure 3 we also display an extended view of NPM-ATI data 
together with the gamm a-ray light curv^j] observed by INTEGRAL satellite around 6:48 UT 



(IMereghetti et al.ll2009l ). Due to the high sensitivity of the SAVNET facility, we can clearly 
see rapid amplitude and phase changes at the corresponding times of the short repeated 
bursts from SGR J1550— 5418. Therefore, we can robustly claim that the rapid changes 
were caused by the short gamma-ray bursts from the magnetar. We listed in Table 1 the 
properties of the SGR short bursts detected by the NPM-ATI VLF propagation path. Al- 
though these magnetar bursts were detected at other SAVNET receiving stations, like PAL 
(Palmas, TO, Brazil), SMS (Sao Martinho da Serra, RS, Brazil), PIU (Piura, Peru) and 
EACF (Antarctica), in this short paper we concentrate on the records from ATI receiving 
station. A detailed comparison of VLF phases and amplitudes observed by the other re- 
ceivers is out of the scope of this paper, and will be reported in a subsequent forthcoming 
article. 



i 



http://www.isdc.iinige.ch/integr al/ibas/cgi-bin/ibas^ics_web.cgi 
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3. Analysis 

To investigate the influence of gamma-ray irradiation in the lower ionosphere, we need 
to characterize the VLF amplitude and phase data. To do this, we estimated the maximum 
amplitude (AA) and phase variations (A0) as well as the recovery timescales of the amplitude 
as follows. First, we extracted the data before and after each burst, and chose a proper 
functional form to represent the baseline level. Most of the baseline levels can be well fitted 
by first-order polynomial functions. When the baseline levels showed curvature, we used 
second-order polynomial functions to represent them. After subtracting the trends from the 
amplitude and phase data, we obtained A A and A0, which are tabulated in Table 1. To 
estimate the typical errors of AA and A0, we made histograms of the residuals, which were 
distributed around with a Gaussian-like form. Therefore, we fitted the histogram using a 
Gaussian, and took the variance as a typical error. 

To quantify the recovery timescales of the VLF amplitude data, we have used the 
function 

Fq 

f(t) = (Baseline level) - — — — — — (1) 

{exp((t - t)/tfaii) + exp((t - to)/t ICV )\ 

where F is a typical amplitude decrease, tf a u is a falling time, and t rcv is a recovery timescale. 
We fitted the data using this function, and determined t rcv for each burst (see Table 1). 



Discussion 



4.1. Amplitude and phase changes 



In Figures 4 (a) and 4 (b), we plot A A and A0 against the gamma-ray fluence from 25 
keV to 2 MeV for each burst (see also Table 1). Although there are not many data points, 
possible correlations between AA and A<f>, and gamma-ray fluences are seen. We note that 
similar correlations were als o reported in the case of X-ray solar flares (IMcRae fc Thomson 
2004! : IPacini fc RauliiJl20~0~6l ) 



We can understand these evidences of correlation in terms of the lowering of the reflection 
height due to gamma-ray ionization. Under a typ ical undisturbed night time condition, VLF 
waves are thought to be reflected at ~85 km (e.g. JCarpenter et al.lll997| ). When gamma-rays 
are directed onto the Earth, they deposit most of their energy in the lower ionosphere, ionize 
the neutral atmosphere there, and produce free electrons. The typical altitude where these 
free electrons ar e produced depen ds on the photon energy. For example, by using Monte 
Carlo simulation llnan et al.l (119991 ) reported that 3 keV and 10 keV photons mainly ionize 
the atmosphere at ~82 km and 60 km, respectively. Similar calculations have shown that 
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gamma-ray illuminatio n increases electron number density below 
photon spectrum (e.g. 



Brownlll973l : lBairdlll974J : lTanaka et al.ll2008h . Consequently, the VLF 



■85 km, depending on the 



radio waves are reflected at a lower altitude than usual, and hence the phase of propagating 
VLF waves is advanced. 

Due to the lack of obser yation, the exact photo n spectrum for each burst was not 
reported so far. But following iMereghetti et al.l (120091 ) we assume the spectral shape of an 
optically thin thermal bremsstrahlung with kT=40 keV. Then, the number of higher-energy 
photons increases as the fluence goes up. Higher-energy photons can penetrate deeper at low 
altitude and increase the electron number density there. As a result, the reflection height 
becomes lower as the gamma-ray fluence increases. 



By treating the propagation of VLF radio waves using the mode theory ( j Wait &: Spies 



1964 ). we estimated the reduction of the reflection height AH from the phase change A(j>. In 
the follo wing, we have assum ed that the lower ionosphere is isotropic and a sharply bo unded 



mediu m (IWait fc Spied Il964j ) . and we have used a phase velocity expres sion given by 



Wait 



( 19591) . Then, the relation between A0 and AH can be expressed as (jlnan fc Carpenter 
1987h 



dAH 



he 



h 
2R f 



+ C„ 



(2n - 1) A 
Ah ' 



(2) 



where d is the length of the disturbed region along the great circle path, / is the wave 
frequency, h is the typical reflection height, c is the speed of light, R e is the Earth's radius, A 
is the wavelength of the VLF radio wave, and n is the order of the waveguide mo de. For a long 
propa gation distance and a normal nighttime reflection height of h ~ 85 km, IWait Spies 
( 119641 ) showed that the second mode (n=2) would be dominant. Therefore, we were able to 
calculate AH from the observed Aip for each burst and these values are tabulated in Table 
1. We note that AH calculated by using Eq. (2) is a rough estimate, and Monte Carlo 
simulations are required to obtain more accurate values. Nonetheless, this gross estimation 
would be meaningful to qualitatively consider the effect of gamma-ray illumination. 

Next, we consider a mechanism for the observed decrease of VLF wave amplitude during 
the gamma-ray illumination. We can understand it on the basis of the altitude dependence 
of the collision frequency u e between electrons and neutral atoms. u e is higher for lower 
altitudes, and it is oft en modeled as u e = 1.816 x 10 11 exp(— 0.15z), where z is the altitude 



measured in km (e.g.. IWait &: Spied 1 1 9641 ) . Consequently, as the reflection height becomes 



lower, VLF radio waves are more attenuated, and hence their amplitude decrease. 

These hints of correlation suggest that it would be possible to deduce a gamma-ray 
fluence from AA and A0. We note that these relations are applicable only for this particular 
VLF frequency (21.4 kHz). There are also another uncertainties which might affect AA 
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and A(f> such as the altitude profile of the ambient electron number density. Nonetheless, 
we claim that the Earth's ionosphere can be used as a new gamma-ray 'detector' and VLF 
data can provide a unique information on incident gamma-ray fluences, even if satellites in 
space were not able to observe it. Therefore, we stress here that this VLF method is a new 
potential technique for monitoring high energy transient phenomena in the universe, once 
we know in advance which source is active. 



4.2. Recovery timescale 



We plot in Figure 4 (c) the recovery timescale of each burst against the gamma-ray 
fluence. We did not find any significant correlation between the two quantities, and all the 
recovery timescales were in the range of 2—12 s. Since the t rcv are longer than the burst 
durations (see Table 1), the observed VLF amplitude and phase time profiles are different 
from the gamma-ray light curves. We also fitted the phase data with a similar function 
of Equation (1) and calculated the t rcv . Again, we did not find any significant correlation 
between gamma-ray fluences and recovery timescales. 

As shown in Figure 2, recovery time profiles can be well represented by an exponential- 
like function of a single parameter t rcv . On the ol 



' Irian et al. 


2007; 


Tanaka et al. 


2008) 



Namely, an initial rapid recovery of a few seconds was followed by a long enduring recovery 
lasting for >1 hour. The long-duration recovery is interpre ted as due to the neutralization 



of positive and negative ions at an altitude below 60 km ( llnan et al.l 120071 ). which means 
that ionization by gamma-rays mainly occurred at such a low altitude. In fact, it is known 
that the photon spectrum and gam ma-ray fluence of giant flar es are much harder and higher 



than those of short repeated bursts (IWoods fc Thompsonll2006l ). Lack of such a long-duration 



recovery in our VLF data suggests that the spectra of short repeated bursts were relatively 
soft compared to that of giant flares. 

As shown above, we interpreted the VLF amplitude and phase changes as due to the 
lowering of the reflection height. In this case, faster recovery timescales are expected for 
larger ga mma-ray fluences , because the electron attachment rate is a negative function of 
altitude (IRowe et al.lll974f ). However, as shown in Figure 4 (c), we did not find such a trend. 
This might be due to d ifferent spectrum f or eac h burst, contrary to what we have assumed 
in this paper following iMereghetti et al.l ( 120091 ). The harder the spectrum, the lower the 
reflection height and the faster the recovery time. Monte Carlo simulations are required to 
confirm this possibility, and the results will be reported in a subsequent article. 
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5. Summary 

We have detected, for the first time, ionospheric disturbances caused by short repeated 
gamma-ray bursts from a magnetar. Amplitude and phase changes of VLF propagating waves 
are correlated with gamma-ray fluences. This can be understood in terms of the lowering 
of the reflection height. While satellites in space cannot continuously observe the whole sky 
due to Earth occultation, the Earth's ionosphere can monitor it without interruption. VLF 
observations provide us with a new method to monitor high energy transient phenomena of 
astrophysical importance. 

Y. T. T. is supported by a JSPS Research Fellowship for Young Scientists. J. P. R. 
thanks FAPESP (Proc. 06/02979-0), CNPq (Proc. 304433/2004-7), and MACKPESQUISA 
funding agencies. 
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Fig. 1. — VLF propagation path from NPM transmitter (Hawaii) to ATI observing station 
(Sao Paulo, Brazil). Also shown are the locations of other four VLF transmitters (NLK, 
NDK, NAA, and NAU). Shaded hemisphere indicates the night side part of the Earth at 
6:48 UT, when the largest burst occurred (see Table 1). The part of the Earth illuminated 
by gamma-rays at 6:48 UT is also drawn by dashed area. 
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5:15 5:20 6:40 6:45 6:50 8:15 8:20 8:25 



Fig. 2. — Amplitude and phase variations of a VLF signal from NPM transmitter (21.4 
kHz), which were observed at ATI (see Figure 1) from 4:00 UT to 10:00 UT on 22 January, 
2009. Lower figures are background-subtracted blown- ups at time ranges during which short 
repeated SGR bursts were detected (see also Table 1). 
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Fig. 3. — (a) Blown-up VLF amplitude data from the NPM-ATI path around 6:48 UT. The 
vertical dashed line shows the time 6:47:57.1 UT, when a relatively large gamma-ray flare 
was observed by INTEGRAL (see also Table 1). (b) Same as (a), but for VLF phase data, 
(c) INTEGRAL/SPl-ACS light curve around 6:48 UT. Note that the peak of the brightest 
burst at 6:48:04.3 UT was probably high er than shown here, due to a saturation problem for 
high count rates ( jMereghetti et al.l 120091 ) . 
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Fig. 4. — (a) Relation between observed ampli tude changes and gamm a-ray fluences (25 keV 
- 2 MeV) measured by INTEGRAL satellite JMereghetti et alJboogh . Values are tabulated 
in Table 1. (b) Same as (a), but for observed phase changes, (c) Same as (a), but for recovery 
timescales estimated by fitting. 



Table 1: Short repeated bursts from SGR J1550-5418 detected by NPM-ATI VLF data 



Start UT 


INTEGRAL 


Duration^ 


Gamma-ray fluence^ 


Amplitude 


Phase 


Recovery timescale 


Lowering of 


on 22 JanJ 


IDt 


[b] 


10~ 5 [erg cm~ 2 ] 


change [dB] 


change [deg] 


of amplitude [s] 


reflection height [km] 


5 17 51.7 


85 


0.45 


1.91±0.01 


-0.33 ±0.03 


2.5±0.3 


7.1±1.2 


8.0±1.0 


5 18 39.5 


93 


1.00 


1.44±0.01 


-0.32 ±0.03 


2.0±0.3 


2.0±0.7 


6.4±1.0 


6 41 02.1 


108 


1.00 


2.35±0.01 


-0.59 ±0.04 


2.2±0.3 


6.3±0.6 


4^5±0.6 
5.8±0.6 


6 44 36.4 


117 


1.75 


1.93±0.01 


-0.91 ±0.05 


2.8±0.3 


5.4±0.6 


6 45 13.9 


121 


1.45 


>4.59±0.01 


-1.8 ±0.05 


6.5±0.3 


5.9±0.3 


13±0.6 


6 47 57.1 


141 


0.35 


1.82±0.01 


-0.74 ± 0.04 


2.2±0.5 


* 


4.5±1.0 


6 48 04.3 


149 


8.15 


>27.76±0.03 


-4.4 ±0.04 


29±0.5 


9.7±0.2 


60±1 


8 17 29.4 


176 


6.20 


6.59±0.02 


-2.4 ±0.1 


9.3±0.3 


12.0±0.4 


14±0.5 



T Taken from 



Mereghetti et al 



(120091 ). The energy range is 25 keV to 2 MeV. * Unable to determine by fitting. 



